function erf(x)
implicit none
real(8) :: erf,x,derf
erf=derf(x)
end function erf
